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Abstract —This paper presents in brief the formulation used for the studies of lateral and lon¬ 
gitudinal dynamics associated with separation and pull out of an ongoing functional stage from the 
spent stage of a multistage launch vehicle using retro rockets as jettisoning system. A fourth-order 
Runge-Kutta method is employed to solve the governing nonlinear differential equations of motion 
to obtain the state vectors of the separating bodies. A statistical method is adapted to account for 
the off nominal design parameters of the bodies undergoing separation and the jettisoning system. 
The tolerances in magnitude of the pull out time, relative velocity, lateral gap and longitudinal clear¬ 
ances between the separating bodies were used to find out the nonoccurrence of collision between the 
separating bodies. This procedure helps in designing and adapting a suitable jettisoning system to 
achieve collision free separation of multistage launch vehicle stages. © 2005 Elsevier Ltd. All rights 
reserved. 
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NOMENCLATURE 


A z 

launch azimuth (deg) 

Re, Rp, Rs 

Earth’s equatorial radius, polar 

h 

altitude of separation (km) 


radius, radius at the surface (km) 

I 

inertia tensor (kg-m 2 ) 

r ,q,p 

body angular velocity along yaw, 
pitch and roll axes (deg/s) 

t 

Time (s) 

X, Y, Z 

body axes yaw, pitch and roll 

mjvf, ms 

mass of the ongoing and spent 
bodies (kg) 

Rext — (Fg ~b 

Ft+ Fj ) 


total external force due to gravity, 
thrust, and jet damping (N) 
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M ex t 

— {Mg -f M t + Mj) 

total external moment due to 
gravity, thrust and jet damping 

h 

distance of thrust location from 
nose tip along longitudinal 
axis(m) 


(N-m) 

d t 

distance of thrust location from 

rj 

position vector in ECI frame (X/, 


vehicle longitudinal axis (m) 


Y 7 , Z 7 ) (m) 


body orientation angles along yaw. 

tog 

position vector of the body CG 


pitch and roll in LPI frame (deg) 


offset (m) 

^GC 

geocentric latitude of the launch 

V 

velocity in the body frame (m/s) 


station (deg) 

m* 

propellant mass flow rate (kg/s) 

© 

launch station longitude (deg) 

T 

T T 

thrust (N) 

distance of thrust location from 

fj- 

Earth’s gravitational constant 
(m 3 /s 2 ) 

ICG 

main body CG (m) 

distance of CG from nose tip 
along longitudinal axis (m) 


Earth’s angular velocity (rad/s) 


1. INTRODUCTION 

Launch vehicles almost always employ multistage using the concept of staging for boosting pay- 
load capabilities. This in turn requires jettisoning the spent stages. The dynamics of separating 
bodies has received the attention of several investigators [1-24]. Chubb [2] has constructed colli¬ 
sion boundaries between two separating stages. Dwork [3] and Wilke [4] provide valuable insight 
into disturbances caused by separation mechanisms in a spinning setup. Considerable amount 
of aerodynamic data including stability derivatives have been generated for space shuttle type of 
configurations involving winged bodies [5,6]. The data so gathered have been utilized for separa¬ 
tion dynamics investigation [8], Waterfall [9] investigated multispring system for separation and 
spinning and nonspinning bodies. Longren [10] analyzed spin stabilized rockets with guide shoes 
and rails constraining the lateral motion. Hurley and Carrie [11] reviewed the genesis of a four 
bar linkage separation system and carried out the analysis for separation of parallel staged shut¬ 
tle vehicles. Subramanyam [13] developed a general model for spring assisted stage separation. 
Biwas [15] investigated several aspects of separation for strap-on from the launch vehicle. The 
aerodynamics of separation is rather complex due to relatively large angles of incidence, intricate 
geometry, and the involved interference phenomenon. Moraes Jr. et al. [16] have developed the 
Navier-stoke’s code, which provided a good agreement with the experiments only on the nose 
part of the ongoing stage. Lochan et al. [17-20] have analyzed the separation dynamics of strap- 
on boosters from the core rocket utilizing the wind tunnel simulation data for dynamic forces. 
Cheng [21] has developed an analytical procedure based on a coupled gas / structure model to 
simulate the fairing separation events. The procedure was validated by comparing the analysis 
results with full-scale payload fairing separation test data of Titan IV launch vehicle. Seongjin 
Choi et al. [22] have developed an efficient three-dimensional aerodynamic-dynamic coupled code 
to simulate the separation dynamics of strap-on boosters in the dense atmosphere. Jeyakumar 
and Biswas [23-25] have provided stage separation system design and dynamic analysis of ISRO 
launch vehicles. 

A typical launch vehicle may involve several separation events such as strap-on separation, 
stage separation, heat shield separation, ullage rocket separation and spacecraft separation (see 
Figure 1). In a multistage rocket configuration the most important event is staging. The staging 
process during the flight of a launch vehicle comprises of a series of events, which takes place 
in transferring the regime of an ongoing stage to that of the next stage. Obviously, the process 
commences from the detection of the burn out of the ongoing stage, and continues with the 
ignition of the next stage and the separation of the spent stage. Ignition of the next stage and 
the separation of the spent stage may takes place concurrently or with a time difference and 
in any order, depending on vehicle configuration and mission requirements, and they may be 
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Figure 1. Separation events of a typical launch vehicle. 


programmed before lift off or on the basis of real time in flight monitoring. During the staging 
process, the orientation control of the vehicle has to be maintained within specified limits - to 
be defined on consideration of controllability, structural integrity and avoidance of recontact. 
Successful staging depends on the proper functioning of the various vehicle systems under the 
environments and disturbances acting on the vehicle. Separation involves pyro-actuated severence 
of the designated vehicle elements and imparting of the required relative velocity between the 
separating bodies through passive or active means. The pyro actuator has to have the capacity 
with specified minimum margin to do the severence and yet it should not generate shock levels 
strong enough to derange the modules of the vehicle. The modules located in the vicinity of 
the severence zone are susceptible to the severence shock. The relative velocity between the 
separating bodies is generated either by releasing a set of precompressed springs or by firing 
auxiliary rockets. In this activity, angular motion is also unavoidably imparted to the separating 
bodies. 

For carrying higher payloads, the launch vehicles are often configured to use solid strapon 
motors. When strapon motors are employed one major concern is the differential thrust among 
the motors. The differential thrust is highest during the tail off period, which is also the zone of 
stage separation. This places high demands on control power plant and hence is to be limited. 
The problems gets complicated when the control power available is also small in nonautonomous 
control systems like the secondary injection or flex nozzle control systems commonly employed in 
solid motors. During the separation event, the vehicle control has to pass from the separated stage 
to the succeeding stage. Where both the stages are solid motors and control is through injection 
of a fluid into the nozzle divergence (secondary injection thrust vector control (SITVC)), some 
special problems arise, during stage separation. In these cases, control is available only during 
the motor operation. Even within motor operation, SITVC effectiveness diminishes for motor 
chamber pressures below a certain value. This situation occurs during the later portion of motor 
tail-off and during the initial phase of ignition. Hence, if the succeeding stage is not ignited before 
the pressure drops below a certain value in the earlier stage, loss of vehicle control occurs. 
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This will result in building up of vehicle orientation errors. This error build-up in the presence 
of atmospheric dynamic pressure loads on the vehicle could prove detrimental to the vehicle 
and the mission. A major factor to be considered in designing the separation system and flight 
sequence is the tail off thrust. Consider the case of separation between two core stages, with lower 
stage being a solid motor. During the lower stage operation, the vehicle control is obtained by 
the stage motor action itself (by SITVC or flex nozzle control). The stage separation should be 
achieved and the upper stage should take over control with out creating large error in orientation. 
The separation should be carried out before the pressure in the booster solid motor (lower stage) 
drops to such a low value that its control is not effective. But if separation is done too early, the 
thrust from the booster sage (tail off thrust) leads to collision between the two separated stages. 
To ensure collision free separation, retro rockets are incorporated on the separated booster. The 
sizing of the retro rockets (thrust levels), booster tail off thrust, dynamic pressure on vehicle, 
vehicle configuration at interstage, separation system force should all be analyzed interactively 
in the separation system design. When the solid motor is followed by liquid engine there is a 
requirement of certain level of positive acceleration before the liquid engine ignition. The ullage 
motor is provided to ensure this level at separation. The ullage motor is ignited before the thrust 
from the lower stage solid motor drops below a certain level. 

This paper presents a mathematical formulation for the analysis of stage separation using retro 
motors as jettisoning systems in a multistage launch vehicle. The problem is modeled under the 
influence of forces and moments of the individual bodies undergoing separation, treating the 
bodies as rigid. A rigid body has six degrees of freedom (three displacement of a point fixed 
in the body and three orientation angles). Each body undergoing separation has twelve state 
variables (namely three displacement of point fixed in the body, three orientation angles, three 
components of velocity vector, and three components of angular velocity vector). Two basic 
types of frames of reference are used in the formulation apart from the usual geo-centric or Earth 
centered inertial (ECI) frame of reference (to compute gravity effect) and the topo-centric or 
launch point inertial (LPI) frame. The first is a local inertial frame (Li-frame) in which the 
dynamics is described and is essentially the body frame of the combined body rotationally frozen 
at the instant of separation command but free to travel at a velocity equal to the combined 
vehicle velocity at separation. This enhances the computation accuracy and makes the results 
of the analysis easy to comprehend. The second type of frame is the body coordinate system 
(B-frame), separately for each of the body undergoing separation. The transformation from body 
frame to local inertial frame and vice versa can be achieved through a transformation matrix in 
which a prefixed sequence of rotation of the Euler angles is used. 

In principle, equations of motion can be expressed in any coordinate system. It is advanta¬ 
geous to express them in the body coordinate system, as the moment of inertia, engine tail-off 
thrust and aerodynamic computations are very simple in the body coordinate system while grav¬ 
ity computation is simpler in geo-centric inertial coordinate system. Transformation from one 
coordinate frame to another is made in the present study as the forces and moments are de¬ 
fined in different coordinate systems. The nonlinear ordinary differential equations of motion 
are solved for a typical retro-rocket assisted separation system using a fourth-order Runge-Kutta 
integration scheme. For separation dynamic analysis, it is necessary to simulate the trajectories 
of the separating bodies taking into account of the forces and moments acting on the bodies. 
To understand the clearances between the separating bodies, critical points will be identified on 
these bodies. The relative distance between the critical points, relative velocity, body angular 
velocity and body orientation angles will be estimated for the identification of collision hazards 
during separation. The trace of the points together with body orientation angles will confirm 
collision free clean separation and build up of any tip of errors on the ongoing functional stage. 
Graphical simulation of separating bodies will show the occurrence of collision (if any) during 
separation. 
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2. FORMULATION 

The physical problem is to simulate the dynamic characteristics of the separating rigid bodies 
under the influence of forces and moments acting on them. Each individual body has twelve state 
variables (namely three displacement of point fixed in the body, three orientation angles, three 
components of velocity vector, and three components of angular velocity vector). In principle, 
equations of motion can be expressed in any coordinate system. It is preferable to express them 
in the body coordinate system, as the moment of inertia, engine tail off thrust, and aerodynamic 
computations are very simple in the body coordinate system. Figures 2 and 3 show the different 
coordinate frames used in the present study. The X/, Y/, and Z/ axes as shown in the Figure 2 
lies in the equatorial plane of the earth. Gravitational force terms are presented in this coordi¬ 
nate system. The Xr, Y t, and Zt axes define the Topocentric or Launch point Inertial (LPI) 
coordinate system. The ECI and LPI coordinates are related through the launch azimuth (A z ), 
launch site geo-centric latitude (3>gc) and launch site longitude (0). The vehicle orientation at 
any time is defined with the Euler angles, viz., yaw orientation angle (ip), pitch orientation angle 
(9), and roll orientation angle (tp) with respect to LPI frame. The body frame coincides with 
LPI frame at the time of launch. It has its origin at the center of gravity (CG) of the individual 
bodies, namely main (ongoing) body and separating (spent) body as shown in the Figure 3. The 
Xm, Ym, and Zm axes with Mcg as origin define the main body coordinate system. The X$, 
Ys, and Z$ axes with 5c g as origin define the separating body coordinate system. The Xl, Yi, 
and Zl axes coincide with body axes of the combined body just before separation. This frame 
moves with a constant velocity same as that of the combined body at the instant of separation. 



Figure 2. Earth centered inertial (ECI) and launch point inertial (LPI) frames. 
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Figure 3. Body coordinate systems. 


2.1. Equations of Motion 


The translational and rotational equations of motion for each body in the respective body 
frame are [7,23-29] 


m 





ext? 


(i) 


—r- + (d X L — Mext- (2) 

at 

Here m is the mass of the body; v = v x i b + v y jb + v z k b , is the velocity; uj — rib + qjb + pkb, 
is the angular velocity, which has components r, q, and p along the body axes, viz. , yaw, pitch 
and roll respectively; t- is the time; F ex t = F x ib + F y jb + F z kb , is the total external force vector; 
L = L x %b + Lyjb + L z k b , is the angular momentum vector; M ext = M x i b + M y j b + M z k b , is the 
total external moment vector. ib,jb, k b are the unit vectors with respect to the center of mass of 
the respective body frame. 
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The components of angular momentum vector are obtained from 



(3) 


Here I xx , I yy , and I zz are the principal moment of inertia about yaw, pitch and roll (X, Y, Z) 
axes, respectively. I xy , I xz , I yx , I yz , I zx , and I zy are the products of inertia with respect to the 
indicated planes, e.g., xy, xz, yz. The symmetry in the products of inertia has been noted. It 
may be noted that the moments of inertia may only be positive, but the products of inertia may 
be either positive or negative. The 3 x 3-inertia matrix on the right of equation (3) is called the 
inertia dyadic [7,10]. 

Equations (1) and (2) can be written in the form 




pv y -qv z ) ( F x 

qv x - rv y ) { F z 

pL y qLz d - Yfj. "I 
rL z - pL x + M y > . 
qL x - rL y + M z ) 



(4) 

(5) 


The kinematic differential equations, which provide the relationship between Euler angles and 
angular velocity components, axe 



cos <p 

— sin >p 

0 " 

i _ _ 

m 

sin <p 

cos cp 

0 


r 

M = 

cos ip 

cos ip 


q 

WJ 

sini p 
simp--- 

, sinip 
cos <p -- 

1 


.p. 


cos ip cos ip 


( 6 ) 


The motion of bodies in yaw plane is important in separation dynamics. Hence, the sequence 
(ip,9,ip) corresponding to the yaw, pitch, roll is followed in the present study, whereas the sequence 
for simulation of flight trajectories adopted is (9,ip, ip). There are singularities in two of the rates 
of the Euler angles 9 and p, which occur in both pitch and roll when ip = 90 deg. These 
singularities in equation (6) are eliminated using L’Hospital's rule [27, pp. 160-163] and one can 
obtain expressions that are finite when ip = 90 deg. In other ways one can resort to Euler four 
parameters approach and enforce a constraint on the parameters with Lagrange multipliers. For 
vertically launched vehicles the yaw angles (ip) will be well below 90 deg. 

Equation (4) to (6) are the governing equations of motion for the individual bodies undergoing 
separation. These are nonlinear ordinary differential equations, which are to be solved numeri¬ 
cally after specifying the components of total external force and moment vectors. The dynamic 
characteristics of the separating bodies are examined here under the influence of the gravity and 
the spring forces and moments. 


2.2. ECI to LPI Transformation 

The transformation matrix [Tel], which transfers the components of a vector from ECI frame 
to LPI frame, is [25-27] 


[Tel] = [e,j]3x3- 


CO 



856 


D. JEYAKUMAR et al. 



Figure 4. Oblate earth geometry. 


Where 

en = — cos A z sin 3 >gc cos 0 — sin A z sin ©; ei 2 = — cos A z sin $gc sin 0 + sin A z cos 0; 

ei 3 = cos A z cos $>gc; 

e 21 = — sin Az sin 3 >gc cos 0 + cos A z sin 0; e 22 = — sin A z sin $gc sin 0 — cos A z cos 0; 

e 2 3 = sin j4 2 cos $gc; 


e 31 = cos $gc cos ©; e 32 = cos $gc sin 0; e 33 — sin ^gc; 


The geocentric latitude of the vehicle position is, $gc — sin -1 (Zj/rj); the geodetic latitude 
of the vehicle position is, $cd = tan _1 {(i?s/i?p) 2 tan$Gc}; Rp is the polar radius of the 
oblate earth geometry; the inertial longitude of the vehicle position is ©j = tan -1 (Y//X/); 
the relative longitude with respect to earth is, ©p = ©/ — flpt; and the earth rotation rate, 
= 7.29211 x 10 -5 rad/s. The altitude of t he vehicle in Figure 4 is h = rj -Rs', the radius of the 

earth corresponding to <3 >gc is Rs = Re /yl + {C R-e/Rp) 2 — 1} sin 2 Tec; the equatorial radius 
R e = 6378.1656 km; and the polar radius, R P = 6356.7838 km. The launch station longitude 
of 0) 80.23621 deg, geodetic launch latitude of (T GD ) 13.73204 deg and launch azimuth of (A z ) 
102 deg E are used in the current studies. 

2.3. LPI to Body Transformation 

The transformation matrix [Tlb], which transfers the components of a vector from LPI frame 
to body frame, is [25-27] 


[Tlb] — [^u'] 3 x3 > 


(8) 
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Where 


&n = cos p o cos 6 q + sin p o sin ipo sin do', &12 = sin po cos ipo; 
613 = — cos po sin 9q + sin p 0 sin po cos d 0 ; 


621 = — sin p 0 cos Oo + cos po sin po sin &o ; ^22 — cos Po cos Po ; 

623 = sin p 0 sin d 0 + cos p 0 sin po cos do ; 

b 3 i = cos ip 0 sin 9 0 ; b 32 = - sin p 0 ', 

833 = cos ipo cos 6 q . 

Here -ipo, do, <Po are the initial Euler angles. It should be noted that the transformation matrix, 
Tlb, which transfers the components of a vector from the local inertial frame (LI) to body frame, 
is nothing but the transformation matrix Tlb having the elements as a function of the sequence 
of rotation with Euler angles ip, 9, and <p. 

2.4. Gravity Force and Moment Computation 

The gravitational acceleration along the ECI frame is computed as [27,28] 

<7EC = 9xie r x + 9Yl^y + 9zie z - -grad U. (9) 


Where the gravitational potential, 


U=- 


9 


r/ 


1 - i hr? (3C 2 - 1) - \ J 3V 3 (5C 3 - 3C) - \j^ (35C 4 - 30C 2 + 3) 


( 10 ) 


V = Re/T, C = Zi/ri\ f / = X[(T x + Yie'y + Zje z , is the vehicle position vector with respect 
to the ECI frame; rj = is the radial distance of the vehicle from the earth center; Re 

is the equatorial radius of the oblate earth geometry (see Figure 4); e x , e y , e z are the unit 
vectors along the axes of the ECI frame; The gravitational constant, /x = 3.9860253 x 10 14 m 3 /s 2 ; 
gravitational harmonics J 2 — 1.0827 x 10 -3 ; J 3 = —2.3 x 10 -6 ; J 4 = — 1.8 x 10 -6 . Since the 
gravitational acceleration (9) is computed in ECI frame, it has to be transformed to the LPI 
frame. The gravitational acceleration along the body frame, g — gxh + 9Y]b + gz^b computed 
using the transformation matrices [Tel] and [Tlb]- The components of the gravity forces in the 
body frame are 


f 9x 1 

1 f 3X1 ) 

J 1 

\ — [Tlb] [Tel] \ 9yi / 

l 9z J 

! { 9zi J 


( 11 ) 


The gravity force and moment are computed as 


F g = mg, (12) 

Mg — fcG X Fg. (13) 


Where fee = XcgH + Ycajb + Zcckb, is the position vector of the vehicle center of mass. In 
general, the center of gravity does not coincide with the center of mass of the vehicle. The 
gravitational forces will cause moment about the center of mass of the vehicle. However, since 
vehicle dimensions are very small compared to the distance from the earth center, the moment 
due to gravitational forces is small compared to other moments. 
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2.5. Retro Rocket Thrust and Moment Computation 

The motion of a spent stage separating from its functional rocket stage by the force of retro 
thrust is considered. The separation impulse is provided by release of the energy stored in the 
relatively low thrust short burn duration solid motor called retro motors. Its location and support 
hardware design is dictated by the geometry requirements, as the jet impingement from the retro 
motor should not affect the highly sensitive payload interface. For each of these thrusts following 
information are required. Time measured from initiation of separation versus thrust in tabular 
form. Thrust action point with respect to CG of separating body on which thrust is acting. 
Mounting azimuth angle with respect to vehicle axis. Nozzle exit area, if separation is in the 
atmosphere and the thrust given is vacuum thrust then the correction is made as follows: [29] 
magnitude of actual thrust = vacuum thrust — (atmospheric pressure * nozzle exit area). 

To compute the actual thrust from a propulsion unit in body frame, a transformation matrix 
is generated between the thrust frame (coordinate frame assumed at the nozzle) and the body 
frame is generated as shown in Figure 5). This matrix is defined by the following sequence of 
anti clockwise rotation of the body frame: about Z-axis through p (measured from Y-axis) the 
azimuth angle defining the propulsion unit location; about X-axis through A- nozzle cant angle 
as shown in Figure 6. 

The thrust vector in the body frame is: 





(A) (B) (C) 


Figure 5. Transformation of retro-rocket thrust frame to body frame. 
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Figure 6. Thrust acting in a retro rocket with respect to main body. 


Propulsion unit thrust location is given by: 

= d? sin pi + dx cos pj + (Icg — It) k- 
Moment due to thrust is given by 


Mt — tt x Ft- 


(15) 


(16) 


2.6. Jet Damping Force and Moments 

Considering a rotating vehicle having an initial pitch rate, and with its motor burning, reveals 
that the mass of gas ejected from the nozzle will have angular momentum. Since the mass is 
isolated from the vehicle once it leaves the nozzle, it is clear that the total angular momentum 
of the vehicle itself must decrease. Consequently, angular momentum is no longer conserved 
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within the vehicle and due to this fact alone; the gyroscopic motion of the vehicle is lessened (jet 
damping effect). The jet damping force and moment axe given by [7,29] 

Fj = — muj x r* 7 ’, (17) 

Mj = tt x Fj, (18) 

where —m is the mass flow rate per unit area (considered a nagative quantity since it decreases 
the vehicle mass). These jet damping forces and moments are computed in the body coordinate 
system. 

2.7. Net Force and Moments 

Using equations (12) to (18), the total external force and moment vectors in equations of 
motion (1) and (2) are written as 


Fext = F g + F t + Fj, (19) 

M ex t = Mg + M T + Mj. ( 20 ) 

With these specified forces and moments, the solution of equation (4) to (6) will give the dynamic 
characteristics of the individual bodies undergoing separation. 

2.8. Identification of Collision between Separating Bodies 

For separation dynamic analysis, it is necessary to simulate the trajectories of the separating 
bodies taking into account of the forces and moments acting on the bodies. To understand the 
clearances between the separating bodies, critical points will be identified on these bodies. The 
relative distance between the critical points, relative velocity, body angular velocity and body 
orientation angles will be estimated for the identification of collision hazards during separation. 
The trace of the points together with body orientation angles will confirm collision free clean 
separation and build up of any tip of errors on the ongoing functional stage. Graphical simulation 
of separating bodies will show the occurrence of collision (if any) during separation [23-25]. 

After determining the velocities of the individual bodies at a time ‘t’, the relative velocity 
between the separating bodies in the body frame is calculated. Applying the inverse transforma¬ 
tion to the relative velocity and integrating with respect to time, the actual separation distance 
between the bodies is found. The integration of the equations of motion is continued till the retro 
rocket is thrusting. To account the variations in the dynamic parameters, a statistical analysis 
is adopted. 


3. STATISTICAL VARIATION OF A FUNCTION 


In launch vehicle separation system design process, it is with the combination of many variable 
parameters. For this purpose a multivariate concept, which involves more variable parameters, 
will be more useful [30]. 

Let the variable function be denoted as T. The variable T be functionally related to l n random 
variables, namely, xi, x 2 ,..., x n . This function is assumed to be continuous and differentiable. It 
can be expressed into a Taylor series about a particular point say xi*,x 2 *,. - .,x n *. Neglecting 
terms of higher order and assuming that the random variables are independent with means at 
xi*,x 2 *, ■ ■ ■ respectively, the expected (average) value of the function T is approximately 


given by 


E(9)n*(xl,xZ,x* 3 ,...,x. 


■;) + ^E Var (*o 


d 2 \E- 


( 21 ) 
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and the variance of the function, T is approximately given by 

Var(¥)«£Vax(*o(|^) 2 . (22) 

The partial derivatives of 4/ in equations (18) and (19) are taken at the point xq*, xq*. . .. ,x n *, 
and thus, are constants. Further more, if the aq are normally distributed, k is approximately 
normally distributed [30], 

If the tolerances of the individual variables all fall with in their specification limit, then the 
tolerances for the function, T with 99.9% confidence level will not exceed E ('I') ± .3^/var (k). 
The square root of the variance of a quantity is the standard deviation of that quantity. 

In the present study, there are 14 variable functions namely velocity components (v x ,v y ,v z ), 
position components ( x,y,z ), angular velocity (r, q.p), body orientation components ( i/),0,4 i>), 
pull out time and lateral gap available at the critical point (nozzle end) during pull out. These 
variable functions are related to the dynamic parameters of the separating bodies like mass, 
moment of inertia, CG offset and separation rocket parameters like cant angle, thrust, ignition 
delay, etc. For the variations in the specified dynamic parameters, the above described statistical 
analysis will be useful to find the expected value and variance of the 14 variable functions for 
the ongoing and spent body. With these values, body rate for ongoing and spent body, relative 
velocity and relative distance between the separating bodies can be estimated. Finally it is 
possible to estimate the tolerances of the physical quantities. 

For the present problem, the random variables xq, x 2 , • • •, x n . are considered as the specified 
dynamic parameters. The variable functions 4b (i = 1 to 14) correspond to velocity components 
(v x ,v y ,v z ), position components ( x,y,z ), angular velocity components ( r,q,p ), and the body 
orientation angles (i/;, 0, </>), pull out time and lateral gap available at the critical point (nozzle 
end) during pull out. The values of v Ib are obtained for the average values of dynamic parameters. 
To find the values of first and second derivatives of 4b with respect to the dynamic parameters 
the following procedure is adopted. 

For each variable function 4b, the functional values generated by specifying a value of dynamic 
parameter close to its average and the rest of the dynamic parameters are specified by their 
average values. The generated data is fitted in to a cubic polynomial of the varied dynamic 
parameter. The first and second derivative of the function with respect to that dynamic parameter 
is found out from the cubic polynomial at its average value. This procedure is repeated for all the 
variable functions and dynamic parameters to obtain the first and second derivative of variable 
functions with respect to each dynamic parameter at its average value. Considering variances in 
the dynamic parameters and the first and second derivatives of variable functions with respect 
to those dynamic parameters equations (21) and (22) give the expected value of the functional 
and its variance. 


4. NUMERICAL RESULTS 

Studies are made on the lateral and longitudinal dynamics associated with separation and pull 
out of an ongoing functional stage from the spent stage of a multistage launch vehicle using 
retro rockets as jettisoning systems. Figure 7 shows the location of retro motors placed in pair’s 
diametrically opposite and the specified pull out length to be cleared to achieve the collision free 
separation. The specified pull out length is 4.6 m. Figure 8 shows the magnitude of actual thrust¬ 
time curve for the retro motors. Tables 1 and 2 gives the trajectory and dynamic parameters of 
reference [23] considered for the separating bodies at the time of separation. Table 3 provides 
the dispersion parameters considered in the studies for the propulsion units. A fourth-order 
Runge-Kutta integration scheme is used to solve the nonlinear ordinary differential equations 
with a fixed time step size of 0.001 s. Table 4 provides the expected value and the variance of the 
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Figure 7. Pull out geometry of the separating bodies and the location of retro rockets 
in a typical launch vehicle. 



71me(s) 

Figure 8. Thrust versus time of retro rockets. 


Table 1. Trajectory parameters at the time of separation. 


Time of separation ‘t’ (s) 

300 

Altitude of separation ‘h’ (km) 

140 

ECI Position Xj, Yj, Zj (m) 

460290.075, 6354251.31, 1432314.34 

Initial body rate r, q , p (deg/s) 

2 for ongoing body and 2.1 for spent body 

Initial orientation \Fo, &o, 4 >0 (deg) 

0.0, 74.5, 0.0 
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Table 2. Dynamic parameters of the separating bodies. 


Parameters 

Ongoing Body 

Spent Body 

Mass (kg) 

34000 ± 100 

12000 ± 100 

Z (mm) 

31000 ± 100 

16000 ± 100 

Y (mm) 

30 ±10 

50 ± 10 

X (mm) 

30 ±10 

50 ± 10 

Izz (kg-m 2 ) 

32000 ± 10% 

50000 ± 10% 

Wy (kg-m 2 ) 

400000 ± 10% 

700000 ± 10% 

Ixx (kg-m 2 ) 

400000 ± 10% 

700000 ± 10% 

Ixy (kg-m 2 ) 

-50 ± 100% 

-150 ± 100% 

Ixz (kg-m 2 ) 

-50 ± 100% 

-150 ±100% 

Iyx (kg-m 2 ) 

-50 ± 100% 

-150 ± 100% 

Iyz (kg-m 2 ) 

-50 ± 100% 

-150 ±100% 

Izy (kg-m 2 ) 

-50 ± 100% 

-150 ± 100% 

Izx (kg-m 2 ) 

-50 ± 100% 

-150 ± 100% 


Table 3. Sensitivities of the propulsion units. 


Parameters 

Dispersion level 

Jettisoning rocket thrust dispersion (%) 

±3.5 

Jettisoning rocket ignition delay (ms) 

100 ± 20 

Jettisoning rocket cant angle (o) 

±0.3 

Differential thrust between the core engines (s) 

±5kN 

Thrust line offset between the core engines (mm) 

±0.3 


Table 4. Variations in 12 degrees of freedom for the separating bodies. 


Physical 

Quantities 

Ongoing Body 

Separated Body 

—3cr Value 

Expected 

Value 

+3(7 Value 

—3 a Value 

Expected 

Value 

±3<x Value 

v x (m/s) 

10.53984 

10.84993 

11.16001 

9.77148 

10.15609 

10.54070 

v y (m/s) 

-2.69797 

-2.62983 

-2.56169 

-2.33249 

-1.92549 

-1.51848 

v z (m/s) 

-4.05700 

-3.93076 

-3.80452 

-13.02184 

-12.88692 

-12.75201 

x (m) 

6.52351 

6.92109 

7.31868 

5.93653 

6.31467 

6.69281 

V ( m ) 

-1.87954 

-1.78350 

-1.68746 

-1.31949 

-1.1757 

-1.03244 

z(m) 

5.25069 

5.39504 

5.53940 

-14.95139 

-14.78274 

-14.61409 

r (deg/s) 

1.41327 

1.41340 

1.41353 

-0.06321 

1.51740 

3.09800 

q (deg/s) 

1.41616 

1.41631 

1.41646 

0.83178 

1.51716 

2.20254 

V (deg/s) 

0.59143 

0.59272 

0.594 

0.55548 

0.60112 

0.64676 

(deg) 

1.70384 

1.75528 

1.80673 

0.80480 

1.88540 

2.96600 

6 (deg) 

1.70256 

1.75409 

1.80561 

1.38802 

1.87811 

2.36820 

<t> (deg) 

0.68658 

0.70813 

0.72969 

0.67446 

0.71442 

0.75439 


12 degrees of freedom of the separated bodies. Table 5 provides the tolerances of the parameters 
of the separating bodies with 99.9% confidence level. It is noted from Table 5 that the expected 
pull out time is 1.24 s and its tolerance is from 1.2 s to 1.276 s. The initial longitudinal distance 
just before pull out between the CG of the separating bodies is 15.561 m. The expected relative 
longitudinal distance between the CG of the individual body at the specified separation time is 
20.19588 m; hence, the actual pull out completed at the specified separation time is 4.63488 m. 

The magnitude of the body rate is: 


\/r 2 + q 2 + p 2 . 


(23) 
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Table 5. Tolerances of the parameters of the separating bodies with 99.9% confidence 
level. 


Parameters 

—3er Value 

Expected Value 

+3<t Value 

Pull out time (s) 

1.20416 

1.24056 

1.27696 

Relative velocity (m/s) 

8.81871 

9.00902 

9.19933 

Relative longitudinal distance (m) 

19.97342 

20.19588 

20.41835 

Lateral gap available (mm) 

646.70596 

801.44444 

956.18293 

Body rate for ongoing body (o/s) 

2.08646 

2.08685 

2.08724 

Body rate for the spent body (o/s) 

1.09493 

2.26818 

3.44127 


| 



Time(s) 


Figure 9. Relative velocity of the separating bodies with time. 


The expected body rate of the ongoing stage is found to be 2.08685 deg/s, whereas in the case 
of the spent stage, it is found to be 2.26818 deg/s. 

The relative velocity (v R ) of the separating bodies is: 

v R = ^ ( v x i - v x2 ) 2 + {Vyi - v y2 ) 2 + ( v z i - V z2 ) 2 . (24) 

The relative distance ( d R ) from the CGs of the separating bodies is: 

d R = \](or - x 2 ) 2 + (yi - Vi ) 2 + (zi - ^ 2 ) 2 . (25) 


Here (v xl ,v yl , v zl ), (v x2 , v y2 , v z2 ), (x^yu z r ) and (x 2 ,y 2 , z 2 ) are the components of the velocity 
vector and position vector of the ongoing and spent body which were transformed to common 
body inertial (CBI) frame. Figures 9 and 10 show the variation of the relative velocity and 
relative longitudinal distance between the separated bodies with separation time. It can be seen 
from Figure 9 that the relative velocity gradually increases upto 8 m/s till the retro action time. 
The average velocity is found to be 4 m/s. It can be seen from Figure 10 that the relative 
longitudinal; distance increases with separation time up to 6 sec due to retro action time and 
decreases with time due to the spent stage tail-off thrust. Hence to avoid collision between the 
separating bodies, the ongoing functional stage should be ignited sufficiently before 10 sec. The 
stage sequencing should be made accordingly. Figure 11 shows the animation of the stage pull 
out at the beginning of separation and Figure 12 shows the animation of the separated bodies 
at the end of separation. In the present analysis nozzle end diameter is 1.6 m. The spent body 
diameter is 4 m; the clearance available (excluding the attachments) for the pull out is 1 m. From 
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Figure 10. Relative longitudinal distance between the separating bodies with time. 



Figure 11. Animation of the stage at the starting of separation. 



Figure 12. Animation of the stage pull out at the end of separation. 


this study, it is noted that the lateral gap consumed during the nozzle pull out of 4.6 m is 200 mm 
due to the variations in the dynamic parameters of the separating bodies and the jettisoning 
systems. It is noted from the results in Table 5 that the expected lateral gap between the nozzle 
end and the burn out body is 800 mm. Hence a collision free separation of the bodies can be 
expected. 


5. CONCLUSIONS 


Dynamic studies are made on the stage separation process and pull out of ongoing stage nozzle 
from the spent stage of a multistage rocket carrier using retro rocket motors. A statistical method 
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is followed to examine the influence of the design variables on the separating bodies. It is noted 
that the lateral gap consumed is about 200 mm against the available gap of 1000 mm during 
the pull out of 4.6 m ongoing stage nozzle from the spent stage, which ensures a collision free 
separation. 
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